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Abstract — We analyze the achievable rate in interference- 
free wireless networks with physical layer fading channels and 
orthogonal multiple access. As a starting point, the point-to-point 
channel is considered. We find the optimal physical and network 
layer rate trade-off which maximizes the achievable overall rate 
for both a fixed rate transmission scheme and an improved 
scheme based on multiple virtual users and superposition coding. 
These initial results are extended to the network setting, where, 
based on a cut-set formulation, the achievable rate at each node 
and its upper bound are derived. We propose a distributed 
optimization algorithm which allows to jointly determine the 
maximum achievable rate, the optimal physical layer rates on 
each network link, and an opportunistic back-pressure-type 
routing strategy on the network layer. This inherently justifies 
the layered architecture in existing wireless networks. Finally, 
we show that the proposed layered optimization approach can 
achieve almost all of the ergodic network capacity in high SNR. 

I. Introduction 

In wireless communication networks, channel impairments 
such as fading, shadowing and path loss limit the capacity. 
Determining the capacity region and corresponding achievable 
strategy for general multiterminal networks is a long-standing 
open problem. An outer bound for the capacity region is 
known to have cut-set interpretation, but this cut-set bound is 
not always achievable. For example, there is no known scheme 
to achieve this outer bound for the simple relay channel. In 
contrast, in wireline networks with orthogonal channels, the 
cut-set bound is achievable even in the multicast scenario, 
where multiple destinations demand the same information 
from several sources. In this case, the capacity of wireline 
networks can be achieved via a separation of physical layer 
channel coding and higher layer network coding. Furthermore, 
in [1] it is shown that the cut-set bound is tight in wireless 
erasure networks. Also, related work in [2] shows that the 
capacity of networks with deterministic channels and broad- 
casting has a cut-set formulation. The optimal rate trade-off 
between physical and network layer is found analytically in [3] 
for a point-to-point channel under an overall delay constraint, 
and by simulation in for minimum energy consumption 
in cellular broadcast networks. However, similar achievability 
results do not yet exist for more general wireless networks 
with fading channels. 

In this paper, we study the achievable rate and optimal 
physical layer rate allocation for a special class of wireless 
networks, where for each node there is no interference be- 
tween the received signals from its neighbors. This could be 
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realized by using orthogonal access strategies such as time-, 
frequency-, or code-division multiple access. Further, the chan- 
nels between each node and one of its neighbors are modeled 
as slow fading channels. For the case where channel state 
information (CSI) is only known at the receiver, each node 
simply transmits at a constant rate regardless of the fading 
state. Therefore, each node's transmitted information gets lost 
when the transmission rate is high so that the current channel 
cannot support it, which results in outage. For the point-to- 
point case we present a throughput-efficient coding scheme 
by partitioning the transmit power across multiple virtual 
users with different rates. The corresponding information is 
then transmitted using superposition coding. This approach 
outperforms the constant rate scheme in terms of achievable 
rate. For an interference-free network with fading channels we 
consider the outage probability as the link erasure probability 
in [1]. Different from [1] where the erasure probability is 
given, the outage probability depends on the transmission rate 
at each node. We first obtain the achievable rate of this network 
and its outer bound based on a cut-set formulation. Further, by 
using a flow formulation, we derive a distributed optimization 
algorithm based on a modified dual decomposition approach 
to maximize the achievable rate. As a byproduct, the algorithm 
also returns the optimal physical layer rate at each node, and 
an optimal opportunistic back-pressure-type routing strategy 
on the network layer Finally, we bound the gap between the 
achievable rate and the outer bound. The result reveals that 
the proposed layered structure can achieve almost all of the 
ergodic network capacity in high SNR, which justifies the 
pervasive layered architecture in existing wireless networks. 

II. Preliminaries 

A. Network Model 

We model the wireless network by a directed acyclic 
graph Q={'V,£), where V is the vertex set and £ is a 
directed edge set. For each node ieV, Mo{i) and Afi{i) 
denote the set of in-neighbors and out-neighbors of i, i.e., 
Miii)={j\{j,i)e£}, Mo{^)^{J\{^,J)e£}. A cut for x,yeV 
is a partition of V into two sets Vx and Vy—V^ such that 
x^Vx and yEVy. Furthermore, we define V* to be the set of 
nodes in Vx that have at least one outgoing edge in V^, i.e., 
V*={i\3j, such that We focus on acyclic graphs in 

the following. 

In this paper, each edge represents a memoryless 

Gaussian channel from node i to node j. Also, we assume 



that each node has a unit bandwidth. Let Xi be the transmitted 
signal by node i and i/ij be the received signal at node j from 
node i, where the average power of Xi is Pi, e.g., = 
Pi. We have the channel model 

yi,j^Vhi,jXj+Vij, (1) 
where Vhi,j is the channel fading between node i and node j 
and Vij is the additive white Gaussian noise at node j for the 
signal reception from node i. Without loss of generality, we 
assume that Vij has zero mean and unit variance. The variance 
crfj of hij is given and normalized such that 0<(Jij<l- The 
channel gain Vhij is assumed to be constant over each packet 
of size n and to vary identically and independently between 
different packets. We consider a Rayleigh fading channel, 
where the probability density function (pdf) of hi j is 



fih^,J) = — e 



(2) 



and hi,j and hi^r are independent Vj,j'GA/'c)(i)j'7^j'. How- 
ever, the results extend to arbitrary fading distributions and 
correlated fading channels in a straightforward way. 

In the following we address multicast problems where a set 
of sinks I?={(ii,...,d|D|}cV demands all of the information 
from a set of sources. In this paper, we focus on the case 
of a single source sGV. The case of multiple sources can be 
generalized in the same way as in [1]. 

Our encoding scheme is given by layered scheme based 
on a product code [5], where a horizontal channel code is 
used for hop-by-hop physical layer error correction and a 
vertical erasure correcting code is used for end-to-end error 
correction on the network layer. This provides an implicit 
block interleaving as an outage for a physical layer block leads 
to individual symbol erasures on the network layer. 

III. Achievable Rate in Point-to-Point Channels 

We first study the point-to-point channel, where the system 
model y=\/hx+v is as described in the previous section but 
with the subscripts i,j removed for brevity. Our objective is 
to maximize the average throughput by optimizing the rate 
allocation for both network and physical layer 

When CSI is only known at the receiver the transmitter 
cannot adapt its transmission rate according to h. In this case 
the capacity is [6] 

C=£h|ilog(l+ftP)|, (3) 

if maximum likelihood decoding is used [6]. Simple encod- 
ing and decoding techniques for AWGN channels cannot be 
applied directly. Without transmitter CSI, one method is to 
let the transmitter transmit at a constant rate R. Therefore, if 
R<C{h) — ^\og{l+hP), the receiver can successfully decode 
the transmitted signal. To maximize the average throughput, 
we need to solve 



maxi? 
R /22n-i 



f{h)dh—maxRe p^'-' —maxF{R), 



(4) 



where we have used the Rayleigh fading pdf in (|2|l. In (|4|l, 
F{R) denotes the overall rate, R the rate on the physical layer. 



with respect to R and setting the result to zero, we obtain the 
optimal physical layer rate as 

where W{-) is the Lambert W function, the inverse function 
of g{w)—we^ . By substituting R* into F{R) in we obtain 
the maximal average overall throughput as 

2^'^* -1 1 - 1 1 



(6) 



We now discuss a coding scheme which is able to improve 
the throughput compared to the constant rate case. The idea 
was first presented in [7] and is based on superposition coding 
by introducing an infinite number of virtual users at the 
transmitter We then propose a modification of this approach 
by considering the practically more relevant case for a finite 
number of virtual users. 

Suppose that each of the users transmits with a power dz on 
the physical layer At the decoder we can employ successive 
interference cancellation to decode the virtual users' messages, 
where the virtual user with a higher interference level z, 
0<z<P, is decoded first by treating all the virtual users below 
as interference. Thus, different virtual users incur different 
interference levels z at the receiver In particular, the virtual 
user at interference level z transmits at rate r{z)dz, where r{z) 
denotes the marginal rate. Note that at channel fading state h, 
a virtual user with power dz and Gaussian interference 1+hz 
can achieve the rate 

1 h 



1 ( hdz \ dz 



-dz. 



(7) 



21n2 1+hz 

Therefore, the message sent by the virtual user at interference 
level z can be decoded if 



2in2 1+hz —^i^)^ ^'^'^ '^he virtual 
users with higher interference levels are decoded successfully, 
if TTT-TT 1 1^ , >riz'), yz'>z, which is denoted as Condition A 

21n2 1+hz' — V /' — ' 

in the following. If 21772 i+hz — ^(-^) has a unique solution for 
z between and P and r{z) is a strictly decreasing function. 
Condition A is satisfied. Then, the probability that the virtual 
user at interference level z can be decoded is given as 

P.(K^),^)=P^(^^>.W). (8) 

Note that in this case the network layer rate is obtained 
by integrating ([8]l over all z, 0<z<P. Thus, the achievable 
average overall throughput for Rayleigh fading is 



Fir{-))-- 



Vd{r{z),z)r{z)dz^ 



Hz) 



r{.).w^ r(z)dz. (9) 



We need to maximize (|9|l over all possible marginal rate 
functions 7-(-) that satisfy Condition A. Clearly, the marginal 



rate function r{z) = 



satisfies Condition A, in 



21n2(:j,^+2:) 

which case (|9]l becomes the single fixed rate problem (|4]i. This 
means that the superposition approach contains the fixed rate 
scheme as a special case. 

By maximizing ^ for each z individually, we obtain the 
optimal marginal rate function as 

^01^0^2 -2 (l+2g-'z-\/l+4a2g^. (10) 



(21n2)c 



and the term exp(— 2— — ji) corresponds to the network layer We can see that 



21n2 1 + hz ' 



(z) has a unique solution 



rate as a function of R. By taking the derivative of F{R) 



and that r*{z) satisfies Condition A. Therefore, r*{z) is 



the optimal solution. By integrating r*{z) over all z, 0<z<P, 
we get the optimal physical layer rate. 

We now propose a modification of the above scheme for 
a finite number of virtual users. For example, for two virtual 
users we choose the marginal rate function as 

— 1 ^, ifO<z<aP, 



r{z)-- 



22H2 -1 



if aP<z<P. 



(11) 



To satisfy Condition A, we must have R2<Ri- Substituting 
(fTTT i into (|9]l, from ( fTOl i we obtain the physical layer rate 



F(i?i,7?2,«) = 7:log 1+q(2 



+ 2 log I 1- 



(1-a) 



(12) 



Fixing a, we can find the optimal i?i(a) and R2{a) by 
maximizing the two summands in ( fT2] i separately. Substituting 
Ri{a) and R2{a) back into (fT2] |. we can find the optimal a 
by a line search. We show experimentally in Section V that the 
performance for two virtual users is not far from the infinite 
user bound. The approach extends to the case with more virtual 
users in a straightforward way. 

IV. Achievable Rate in Wireless Networks 
In this section, we consider general wireless networks 
without interference. As above, we assume that CSI is only 
available at the receiver 

A. Achievable Rate Region and Upper Bound 

For the sake of brevity we only consider an extension of 
the fixed rate case in Section Hill results similar to those 
presented below can also be obtained for the virtual user-based 
scheme. Each node i simply broadcasts to its neighbors at a 
common rate Ri due to the lack of CSI at the transmitter The 
probability that node j£j\fo{i) can receive the packets sent 
by node i is given as 



f{hi^j)dhij=e 



(13) 



where we have used the Rayleigh pdf from Due to 
the ergodicity of the channel, we can associate each network 
layer packet transmitted by node i with an erasure probability 
eij{Ri) — l~Vi.j{Ri) at node j. Therefore, at the network 
layer, we see an erasure network, which is similar to that in 
[1]. It can be readily verified that the results in [1] can be 
extended to our case with only a slight modification in the 
proof, which leads to the following Theorem. 

Theorem 1. The capacity of the wireless network with fixed 
rate transmission and Rayleigh fading on each link is given by 
the minimum value of the cuts between the source node and 
any of the destinations, i.e.. 



C— max min > 



fi{Ri)Ri 



X 1- Yl e^,,{R,)\dR, 



subject to 



r+oo 

/ f^{Ri)dR, 
Jo 



(14) 



In Theorem 1, the function fi{Ri) plays the 
role of time sharing between different rates. Since 

'^nAAc.(i)}'''j(^«)) 1^ generally not a convex 
function in Ri, time sharing may achieve some rate that 
is not achievable without using time sharing. Note that a 
similar result as in Theorem 1 can also be obtained for the 
virtual-user based scheme from Section [III] 

Since Theorem 1 only addresses the fixed rate transmission 
scheme for the receiver-only CSI case, it is instructive to 
also consider the upper bound on the achievable rate for any 
transmission scheme. Define ^j = {Yij\iEj\fj{j)}. By using 
the cut-set bound from [8, Theorem 14.10.1] we obtain 



|X^=VvV,cV,sgV, andOnV:/0, 

(15) 

where X^^ ={X,\i€Vs}, Y^^ ={Y^\jeV^}, i/[v„v|] = 
{Hij\i£Vs:j&V^,{i,j)e£}, and Ht^j denotes the 
corresponding random variable for hij. As CSI is known 
at each receiver, we can interpret the CSI as an additional 
received quantity at the receiver Since we assume that 
different nodes have orthogonal channels and that the network 
is acyclic, an upper bound on the network capacity of receiver 
CSI only can be obtained from ( fTSl l as 



C< mm 

{VslVsCV.sev^ 



ilog l+P. ^'A} (16) 

\ jevjnAAo(i) /J 
Note that the upper bound is difficult to achieve in a general 
network because the receivers in the cut-set are assumed to 
have full cooperation to obtain the cut-set rate. In particular, 
this holds for an achievable scheme. However, in Section [TV-CI 
we will show that the achievable rate is not far away from this 
upper bound in high SNR. 

B. Opportunistic Routing and Distributed Algorithm 

To solve (O directly, we need to consider all the cuts which 
may be exponential in the number of nodes in the network. As 
an alternative, we consider a flow-based formulation where the 
constraints are only polynomial in |£| and |V|. At each node 
i in the network let the rate on each outgoing link be denoted 
for each destination node d^V. Then, in order to find 



^,3 



the maximum multicast rate, we need to solve 



max C 



s.t. 



E 



E 



C, if i—s, 
— C, if i=d, 
0, otherwise, 



E^M< / f^(Ri)R^\l-Y[e^,,{R^)\dR,yZ(ZMo{i 



f,{R^)dR,=l, 



(17) 



where xfj is the flow rate on link for destination dEV 
and the vector x contains the flow rates for all d and for all 
edges in the network. The right hand side of the last constraint 
in ([TtI i represents the total amount of information that can be 
decoded by at least one node in Z. It can be readily verified 
that (O and (fTTI l have the same optimal value. 

Note that ([TtI i is similar to the problem in [9], where a 
dual decomposition based cross layer optimization is proposed. 



Our problem differs from that in [9] in the realization of the 
physical layer rate, which also results in a different routing 
scheme compared to [9] as shown in the following. By 
adopting a similar dual decomposition approach as in [9], we 
can write the Lagrange dual of ( fTTl ) as 



o(q) 



=max 

x,/,C 



d 



S.t. 



E 

4-00 



f^{R^)R^ \ 1 



L 



MR,)dR,=l, 



dRiyZ<zMo{i 



(18) 



where xf is defined as the right hand side of the first constraint 
in dTTI l. and qf denotes the Lagrangian multiplier at node i for 
destination d. The collection of all qf is given by the vector 
q. Due to the convexity of ( [TtT i there is no duality gap. 

To find £'(q) in ( fTsT l. we find that the maximization problem 
separates in two individual maximizations over x and /, resp., 
and C for given q. In order to update C we use a primal 
subgradient method. Let C* be the value of C at the t-th 
iteration. We update C using 



C'+7 1 



-E^ 

dec 



(19) 



where [•]+ denotes the mapping to non-negative numbers and 
7*>0 is a stepsize. For the second maximization over x, let 



2J 



tt'^ be a permutation of jGAfo{i) such that 
satisfies Wi ^d>Wi ^d>-- ■>w.i _d . Due to the polymatroid 

structure of the constraint in ( flSl l. ( flSl l is maximized when 

MR^)R^(l-e^^,diR,))dR,, 

(20) 



n + oo " i 

-- / MR^)R^ll'^,.f(^^){^~^.^^i{R^ 

Jo 



dR„ 



fc=2,...,!AAc,(j)|. 

We can see that the solution ( |20| | corresponds to prioritizing 
the neighbors of node i for destination d according to wfj. 
Thus, if a higher priority neighbor has received the message 
from node i for destination d, all the lower priority neighbors 
will drop the same message even though they also have 
received it. We can interpret qf as the queue length of the 
flow for destination d at node i and wf^ as the aggregate 
queue length difference between i and j for destination d, 
respectively. Different from opportunistic routing [10] where 
neighbors are prioritized heuristically, we optimize the pri- 
oritization by taking into account the queue information at 
each node. Note that each node's queue length is affected 
by both the incoming traffic rate and the depletion rate. 
The depletion rate implicitly counts the optimization over the 
network beyond node i. Thus, this prioritization approach can 
be seen as special variant of a back-pressure-based routing 
scheme [11]. 

Substituting ( l20l ) into ( fTsl l. we obtain an optimization prob- 
lem which depends only on /i(-). The solution is fi{Ri) = 
d{Ri — R*), where 5{-) is the Dirac delta function, leading to 

|JVo(»)l / 3 

3 = 1 dGT} \ fc = l 



As (l2Tl i is not convex, we can solve it using a line search. 
Note that at different iterations t the chosen rate R* in ( |2T] ) 
can be different due to oscillations of q. Time sharing is then 
realized by averaging R* over all iterations. 

After solving ( fTsl l. we obtain the dual function -D(q), which 
needs to be minimized. Instead of minimizing D{q) directly, 
we use the dual subgradient method to update q. Let q, 
the dual variable at the <-th iteration. We update qf'* 

['li'^-V\T,{j\jeAfom^i,3~^{j\3SAfj:(i)}^ii~ 



be 



9i 



d,t+l_ 



'^i,3~ 



■E 



{j|jeArx(i)}- 



^.d 



Via 

-C 

if i=s, 
if i=d, 
+ 



3,1 J 

Otherwise. 

(22) 

where rf>Q is a constant stepsize, and C* and xf ^ are 
obtained from ( fT9] l and ( l20l l. respectively. The updated qf'*^^ 
is then used in (fTsT l again for the next iteration. By the subgra- 
dient optimization theory and the convexity of the problem, 
the proposed primal-dual subgradient algorithm in (fT9t- (l2T]) 
converges to the optimal solution in (fTsT l when stepsizes 7* and 
77* go to zero as t^oo and converges to a small neighborhood 
of the optimal solution with fixed stepsizes [9]. 

An important benefit of the presented algorithm is that ( fT9l ), 
(I20I 1. and i2T[ can all be solved locally at each node i given the 
queue length information of its neighbors, which only involves 
local information exchange. This algorithm can not only be 
used to determine the maximal multicast rate and the corre- 
sponding physical layer rate on each link but also suggests 
the optimal routing protocol to achieve this rate. We can see 
that a layered structure can be used, where channel coding 
is applied at the physical layer and opportunistic routing and 
erasure correction coding are applied at the networking layer. 

C. Asymptotic Behavior 

In [12] we prove the following theorem which quantifies 
the capacity gap between the achievable scheme proposed in 
Section ITV-BI and the upper bound on the achievable rate from 
Section IRTAI 

Theorem 2. Let C'' {P)=C'' {[Pi,P2,...,P\v\]) be the capacity 
region of the fixed size wireless network without interference 
and C(P) be the achievable rate with fixed rate transmission 



at each node. Further, let <j= 



=maximaxj^j/ afj / af^. Then, 
lim Ct(P)-C(P)<i|O]|V|log(]V|<;)+0.7588|D||V| 



|01^1nln(P. 



(23) 



where tlie limit is taken for each entry of P individually. 
Moreover, limp^ooC^(P)/C(P) = l. 

Theorem 2 indicates that the fixed rate transmission scheme 
can achieve almost all of the ergodic network capacity asymp- 
totically in high SNR. Hence, in high SNR the proposed 
layered structure is close to optimal which justifies the layered 
structure in existing wireless network protocols. However, 
different from existing protocols, in the proposed scheme the 
queue length information determines not only the opportunistic 
routing but also the physical layer data transmission rate at 
each node. 
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Fig. 1 . Comparison of achievable rates for different scliemes in point-to-point 
Rayleigii fading channels. 

V. Simulation Results 
We assume in the following that all the nodes have the same 
transmit power P and that the channel gains have the same 
unit variance. 

Fig. [U compares the achievable rate of different schemes for 
a point-to-point Rayleigh fading channel. The fixed (single) 
rate transmission without transmitter CSI (|5]l is denoted as 
"CSIR, One Rate", the scheme with two virtual users and 
superposition coding in (fTTT l is denoted as "CSIR, Two Rates", 
'CSIR, Infinite Rates" represents the approach for infinitely 
many users and the achievable rate (fTOl i. and "CSIR, Capacity" 
denotes the capacity with receiver-only CSI in (O. Also, the 
capacity with both transmitter and receiver CSI obtained by 
waterfilling [6] is shown as "CSIRT, Waterfilling". We can see 
that the CSIR and CSIRT capacities are the same in high SNR, 
indicating that transmitter CSI is not necessary in this case. We 
also observe that the rate loss of single rate fixed transmission 
to the capacity increases as SNR increases. For two virtual 
users, the achievable rate improves at high SNR, however, still 
with an increasing rate loss to capacity. In contrast, the infinite 
virtual user case exhibits a constant rate loss to capacity. 

We now consider a diamond network with a single source 
node s, a single destination node d, and two relay nodes ri and 
r2- There are no direct links between source and destination 
node and between the relay nodes, resp., and source and relays 
transmit in orthogonal channels. Fig. |2] displays the achievable 
rate and the upper bound (fT6b for different P. As the network 
has only four cuts, we can compute the achievable rate for 
fixed rate transmission by using the cut-set formulation (fl4] | 
directly. We also include the achievable rate using the flow- 
based formulation (T% optimized via the distributed algorithm 
in Section IIV-BI We can see that the flow-based formulation 
achieves the same rate as the cut-set formulation. The rate gap 
between the achievable rate and the upper bound increases as 
P increases but diminishes relatively to the total rate. 
VI. Conclusion 

We have analyzed the achievable average rate in 
interference-free networks along with the optimal rate allo- 
cation on network and physical layers. For the point-to-point 
fading channel, we presented a superposition coding approach 
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Fig. 2. Comparison of rates for a diamondnetworkby maximizing )14t d irec tly 
and by using the proposed distributed algorithm and the upper bound )16t . 

based on multiple virtual users to improve the throughput of 
fixed rate transmission when only receiver CSI is available. For 
networks, a distributed cross-layer optimization algorithm was 
proposed to obtain the maximum achievable rate using a fixed 
rate transmission scheme on the physical layer It was shown 
that the proposed layered approach asymptotically achieves 
almost all of the ergodic network capacity in high SNR. 
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